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Recent spin-Seebeck experiments on thin ferromagnetic films apply a temperature difference AT X 
along the length x and measure a (transverse) voltage difference AV y along the width y. The 
connection between these effects is complex, involving: (1) thermal equilibration between sample 
and substrate; (2) spin currents along the height (or thickness) z; and (3) the measured voltage 
difference. The present work studies in detail the hrst of these steps, and outlines the other two 
steps. Thermal equilibration processes between the magnons and phonons in the sample, as well as 
between the sample and the substrate leads to two surface modes, with surface lengths A, to provide 
for thermal equilibration. Increasing the coupling between the two modes increases the longer mode 
length and decreases the shorter mode length. The applied thermal gradient along x leads to a 
thermal gradient along z that varies as sinh (x/X), which can in turn produce fluxes of the carriers 
of up- and down- spins along z, and gradients of their associated magnetoelectrochemical potentials 
p.f,l, which vary as sinh (x/X). By the inverse spin Hall effect, this spin current along z can produce 
a transverse (along y) voltage difference AV y , which also varies as sinh (xj A). 

PACS numbers: 75.30.-m, 44.10.+i,85.75.-d,85.80.-b 



I. INTRODUCTION 

In principle, a thermal gradient VT can produce a 
spin current P This magnetic analog of the Seebeck ef- 
fect, whereby electric currents are generated by VT, is 
known as the spin-Seebeck effect (SSE). Evidence for the 
spin-Seebeck effect has recently been observed in ferro- 
magnet films with thicknesses d,F ~ 10 nm and lengtsh 
L ~ 10 mm grown on insulating substratesPS When sub- 
jected to a temperature gradient (see Fig. [l^,), a nonzero 
voltage difference AV y across the width of the sample is 
observed; this signal is attributed to an Inverse Spin Hall 
Effect (ISHE) due to an inferred spin-Seebeck-induced 
potential gradient along z. (We employ the magneto- 
electrochemic al po tential p^ t i introduced in Ref. [T| and 
defined in Sec. 
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The magnitude of AV y is observed to 
decay in space over a length much greater than a spin- 
diffusion length. 

The relation between the applied temperature differ- 
ence and the measured voltage difference is complicated; 
the connection is represented by 



AT, 



Equil. SSE „ _ 

> d z T > <9 z /z u 



ISHE 



>AVy, 



(1) 



where AT X is applied and AV y is measured, and "Equil." 
denotes thermal equilibration processes. The present 



work shows the details of AT, 



Equil. 



> d z T, then discusses 



d z T -555> a z/ Lt t 4 and d z ^^ AV y . 

Reference [2| observes the voltage difference AV y along 
y to have a sinh (x/A)-like form along the sample for 
some A = A eX pt, thus indicating a surface effect associ- 
ated with heat input and output. It has been suggestecP 
that this surface effect is governed by magnon-phonon 
thermal equilibration 6 within the sample, which has a 
characteristic length of X mp . However, Ref.[5]argues that 
for permalloy (NigiFeig) this equilibration should yield a 





FIG. 1. The substrate (s, dark gray) and ferromagnetic sam- 
ple (F, light gray) of the spin-Seebeck experiment. Here, (a) 
shows the typical experimental system, and (b) shows the sys- 
tem with a disconnection (scratch) in the sample (but not the 
substrate) of length Id- An external magnetic fiefd B x is ap- 
piied aiong x, and a temperature difference AT X aiong x is 
maintained by a heater and a heat sink. A voltage difference 
AVy across the sampie in the y-direction is measured as a 
function of x by point eiectrodea^ or by Pt wires (not shown) 
deposited on the sampleP^For a scratch length Id = 350 /im, 
Ref. [2] measures a similar signal A V y as for the unscratched 
sample (a). The figures are not to scaie. The heater and heat 
sink, which are piaced at each edge of the substrate afong x, 
are not pictured; see Fig. [2] 
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maximum characteristic length of only X mp = 0.3 mm, 
whereas experiment shows the spin-Seebeck effect to have 
a characteristic length at least an order of magnitude 
larger!^ 

Further, the effect is unchanged for a large 
discontinuity^ (of length l& = 350 fixn) along x (see 
Fig. [TJj) ; both with and without the discontinuity, a sin- 
gle sinh (a;/ A) is measured across the entire length L of 
the system. Clearly the substrate, which is the only phys- 
ical connection between the discontinuous regions of the 
sample, plays an important roleP 

This work studies heat flow due to the excitations re- 
sponsible for thermal conduction - that is, magnons (spin 
waves) and phonons (lattice vibrations) - in this system. 
We employ irreversible thermodynamics to justify and 
extend the fD, two-subsystem approach of Ref. |B]to the 
system of Fig. [T] a 2D (x and z) system with transna- 
tional symmetry along y, which contains three subsys- 
tems: sample phonons (designated by subscript p), sam- 
ple magnons (m), and substrate phonons (s). We con- 
sider that substrate phonons incident on the interface 
directly excite only sample phonons, but not magnons 
(justification for this approximation is discussed below). 
For a non-magnetic sample, the characteristic length of 
the sample-substrate thermal equilibration satisfies (see 
below) 



(2) 



Here, k is a thermal conductivity, d is a length related 
to the thicknesses of the sample and substrate, and Kk 
is the thermal boundary conductance!^ 

In addition to various geometrical lengths, there are 
three different lengths associated with Fig. [TJ the sample 
magnon-phonon equilibration length A mp ; the substrate- 
sample phonon equilibration length A ps ; and an infinite 
length Aoo that leads to the usual linear thermal profile. 
Recall that Ref. [2] observes a sinh (x/ A) profile of the ef- 
fect. If A <C L, then sinh (x/X) can decay too close to the 
boundaries to be experimentally observed. Conversely, 
if A 3> L, then sinh(x/A) will appear to be linear in x, 
which may explain the linear signal observed by Refs. 3 
and |H It is therefore likely that the longer of X ps and 
A mp is the length observed. Moreover, because the re- 
sults are independent of we expect that the longer 
length Ai ^> id and the shorter length A2 <S id- 

When both magnon-phonon equilibration (internal to 
the ferromagnetic sample, and not present for a non- 
magnetic sample) and sample-substrate equilibration 
(not present for a sample with no substrate, as in Ref. [5]) 
are present, the coupling between these two modes fur- 
ther separates their characteristic lengths. That is, the 
longer length Ai and the shorter length A2 are respec- 
tively greater and less than both A ps and A mp . With K m , 
k p , and k s denoting the respective thermal conductivi- 
ties of magnons in the sample, of phonons in the sample, 
and of phonons in the substrate, the coupling is given 
by dimensionless coupling constant that is the product 



of R mp and R ps , which are shown below to be given by 



R 



mp 



1 Rps — 



dpK v + d s K, & 



> (3) 



The thicknesses d s and dp are shown in Fig. [I] Thus, the 
coupling strength (and the increase of the longer equi- 
libration length) is enhanced via R mp if the magnons 
account for an appreciable amount of the thermal con- 
ductivity of the ferromagnet, and is enhanced via R ps if 
the substrate is much thicker or has a much larger heat 
capacity than the ferromagnetic sample. 

Section [TT| employs irreversible thermodynamics to find 
the energy transferred between two systems at different 
temperatures, specifically considering systems that share 
a surface (e.g., the sample and substrate) and systems 
that share a volume (e.g., magnons and phonons in the 
ferromagnet). For heat flow only along x, Section III 
finds the characteristic lengths of the thermal equilibra- 
tion modes, as well as the spatial profiles of the phonon 
and magnon temperatures and heat fluxes. For heat flow 



along both x and z, Sec. IV finds the shape of the spatial 
profile of temperatures and heat fluxes, and numerically 
solves for the characteristic lengths and z-dependence 
of the phonon and magnon heat flux magnitudes. Sec- 
tion [V] compares estimates of the thermal equilibrat ion 
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lengths^ 1 to the observed decay length of AV y . Section 
discusses the connection between the thermal gradients 
found in Sec. |IV| and the magnetoelectrochemical poten- 
tials (which involves the spin-Seebeck effect) and the sub- 
sequent connection to A V V (w hich involves the inverse 
Spin Hall effect). Section |VII| provides a brief summary 
and conclusion. Appendix|A| gives details of the bulk and 
boundary conditions associated with heat flux along both 
x and z, used in the numerical calculations in Sec. |IV| 

It has recently been proposecP^ that electron-phonon 
drag and magnon-phonon drag processes are important 
in explaining the results of Refs. [5HU (The kinetic theory 
of electron-phonon drag is found, for example, in Refs. fTTl 
[12| andfTSl) This work does not consider such effects. 



II. THERMODYNAMICS 

Flow described by thermodynamics is properly given 
by the methods of irreversible thermodynamics. We 
present here a derivation of a result central to Ref. |H1 
which is the basis of Ref. [5J but which is simply written 
in Ref. [H 



A. General Equilibration of Two Systems 

We consider any two systems through which heat and 
entropy (but not matter, quasi-momentum, or momen- 
tum) flow. We later specifically consider energy equi- 
libration between the phonon-magnon subsystems in a 
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ferromagnet (as in Refs. [5] and |H]) , as well as energy equi- 
libration between the respective phonon systems of a fer- 
romagnet and a non-magnetic insulator in contact. 

In two such systems, designated a and (3, the energy 
differentials may be written as 



dE a — T a dS ai dEp — TpdSp, 



(4) 



where T is the temperature and S is the entropy. By 
energy conservation dE a = —dEp, so 



dS a 



dE a 



dS, 







dE a 



(5) 



Since the entropy change must be non-negativep' we 
have 



< S a + Sp 



1 



E n 



Tf3 — Tg 

T a Tp 



For S a + Sp > to hold we must have 
E a — C (Tp — T a ) , 



E a . (6) 



(7) 



where C > 0. That is, by irreversible thermodynamics, 
the energy flux is driven by a difference in intensive ther- 
modynamic quantities. The proportionality coefficient £ 
has units of a specific heat divided by time, and as noted 
below depends either on a boundary conductance (for 
systems that share a common surface) or a relaxation 
time (for systems that share the same volume). 
Specific heats per unit volume (C) are defined via 



£ a — C a T a , ip — CpTp 



(8) 



where e = E/V and V is the volume of the system. Use 
of Eqs. Q and ([£]), and Ep = -E a , yields 



T = 



Tr — T n 



T a = 



T„ — Tp 



(9) 



where r Q = C a V a /C and Tp = CpVp/Q have units of time. 
Then 



AT a « = Tp — T a = - 



Tp — T a 



where we define 



T a p 



TgTp 
T a +Tp' 



(10) 



(11) 



Equation (fO) justifies Eq. (1) of Ref. El 



B. Two Systems Occupying the Same Volume 

Energy conservation in two systems that occupy the 
same volume V (e.g., the phonon and magnon systems 



within a ferromagnet) gives e a = —ip, so that substitu- 
tion of Eqs. ([9| and Q into Eq. M yields 



C a _ Cp _ C 

T Q Tp V 



Then, with Tp = (Cp/C a )T a , equation (ff ) gives 

CgCp 



C a _ Cp 

Tg Tp 



C a + C r i 



13 ^r- 1 



(12) 



(13) 



This is the case studied by Ref. [6l 



C. Two Systems with a Contact Surface 

For two systems in thermal contact over a surface of 
area A (e.g., the ferromagnet and substrate's respective 
phonon systems in Fig.[T|, we write £ = /ix^f^so that 



Eg — —Ep — h K A (Tp — T a ) . 



(14) 



Here \ik is the thermal boundary conductance. Substi- 
tution of Eqs. ([14]) and M into Eq. |8| gives 



(15) 



dgCg dpCp 

h K h K 
where d is the thickness of the material in the direction 



normal to the contact surface. Eq. (If I then gives 



T~a/3 



dgCgdpC'p 



1 



flK \dgCg + dpCp 



(16) 



III. HEAT FLOW IN ID 

We now consider a ferromagnet/substrate system 
where a thermal gradient is applied by a heater at 
x = —L/2 and a heat sink at x = L/2 (see Figure |2j. 
For sample isolation, we take them to be in contact only 
with the substrate. This affects the relative amplitudes 
of temperature and thermal flux in each mode, but does 
not change the mode lengths. 

We now take heat to flow only along the length of the 
materials (the x-direction in Figs. [I] and 2j), i.e., heat 
flow in each system is uniform in the yz-plane (Sect IV 
considers flow along x and z). Conservation of energy, 
with an energy source, is given by 



i + d x j% = 



(17) 



where j £ is the energy (and heat) flux, and S £ represents 
the rate of heat transfer per unit volume from one system 
or subsystem to another. We consider steady state solu- 
tions, so that £ = 0. Further, we take the magnon system 
(m) in the ferromagnet to only transfer energy to/from 
the phonon system (p) in the ferromagnet. Similarly we 
take the substrate (s) to only transfer energy to/from the 
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FIG. 2. An zz-plane cross section of the system under con- 
sideration (see Fig.[l]). The heater and heat sink, represented 
by squares at x < —L/2 and x > L/2, maintain tempera- 
tures Th and Tc, where Th > Tc- For sample isolation, we 
take them to be in contact only with the substrate (s, dark 
gray), and not with the ferromagnetic sample (F, light gray); 
this affects the relative amplitudes of the modes, but not the 
mode lengths. The total heat flux input by the heater at 
x — —L/2 is jo, and a similar heat flux must exit the sub- 
strate at x — L/2. In Section |lll| we further take all heat 
fluxes to be uniform in the yz-plane; this restriction is lifted 
in Section Hvl 



phonon system (p) in the ferromagnet, thereby neglecting 
the magnon-substrate coupling. 

The rate of energy transfer per volume (V = Ad) be- 
tween substrate phonons and sample phonons (an energy 
source S) is found from Eq. ( 14 ) as 



IK 



(T s -T p ), S, 



p— >S 



h K 



(T p -T s ). (18) 



Here S e A ^ B is the volume rate of energy transfer from 
system A to system B. This energy transfer is in the 
form of a source only because here we take the heat flux 
to be only along x. When we include heat flow also along 



z in Sec. IV the substrate-sample phonon energy transfer 
is properly treated as a heat flux along z. 

The volume rate of energy transfer between the 
magnons and phonons in the sample is found by sub- 
stitution of Eqs. Q and (12) into Eq. which gives 



-SI 



(19) 



Here we have used Eq. ( 12 ) to replace C p /t p with C m /r m . 
Applied in turn to the substrate, magnons, and phonons, 
Eq. (17) gives 



d x f x ° = ^(T p -T s ). 



(T m - T p ) , 



d*3l* = ^{T.-T p ) + ^{T m -T 9 ). 



(20) 
(21) 
(22) 



As usual, for each subsystem we take the he at flux to 
be proportional to the gradient of temperature j 1 * 15 *^ so 



(23) 



Here k > 0, i.e., heat flows from hot to cold. We have ne- 



glected cross-terms in Eq. ( 23 ) , where gradients of other 



intensive thermodynamic quantities also cause a flux; we 
discuss these cross-terms in further detail in Sec.lVIl Sub- 



gradient of Eq. ( 23 ) in turn gives 



stitution of Eqs. (20), (21), and (22) into the linearized 



' U; *'%T a =T p -T„ 



-k f) 2 T - 

l% p u x- L P 



h K 

K/ m T n 



On 



ri 2 T — T —T 

C n 



(24) 
(25) 
(26) 



A. Characteristic Lengths 



We denote the inhomogeneous parts of T s , T p , and T m 
with primes. They all vary as e ±qx , so the characteristic 
length is A = q^ 1 . Then, solving Eqs. (24) and (25) for 
T' s and T' m yields 



T 
p 



d s K s 

flK 



tL 



T 
p 



^m^m \ 9 
-FT- ) 1 



(27) 



Substitution of Eq. (27) into Eq. (26) gives 



2 n i< 
-K P q = -r- 
dp 



'-q 



(28) 



This is cubic in q 2 . One solution is q 2 ^ = A^ 2 = 0, 
corresponding to the usual linear temperature profile, for 

which t' s = t; = T' m . 

We define the inverse lengths q mp = A~p and q ps = 
Apj. 1 , the former associated with magnon-phonon equili- 
bration within the ferromagnet and the latter associated 
with substrate-sample phonon equilibration. They sat- 
isfy 



2 C m I K m -\- Kp 

Imp = —~ I ~ ~ 

>m \ r^m^p 



2 _ 
5 Qps — 



I'K 



dpKp + d s K s 



(29) 



They are the inverse lengths of the modes when 
the magnon-phonon system and the substrate-sample 
phonon system do not interact. Then for q 2 ^ 0, equa- 
tion ( 28 ) can be written as 



n — 4 — 2(2 , 2 \ , / 2 2 2 2 n p \ 

u y Q \9mp T q ps ) ' \9mp < Zps qmpqps^mpJ^ps ) j 

(30) 

where the dimensionless ratios R mp and R ps are given by 
Eq. pj). The solutions are 



^2 Qmp qps _j_ If q^np qps ^ > . > 



9(1,2) 



qmpqps^-mpRps i 

(31) 
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where q\ is associated with the minus sign, so that q\ < g 2 
and Ai > A2. 

We now consider two extreme cases. If there is no 
substrate (or if hx — > 0), then 



B. Thermal Profile and Fluxes along x 



We write the phonon temperature in the ferromagnet 



as 



\q\ -> q mp = 




(32) 



which on use of Eq. (13) reproduces the result of Ref. 6 



(which employs A for q). If there is a substrate but no 
magnons (or r m —¥ 00), then 




d s K s + d,FKp 



(33) 



T p = T Q + ax + [ T i sinh ili x ) + T 7 cosh (Qj x )] » 

7=1 

(34) 

where To, Tf , T£ , T x & , and T| are temperatures, and a is 
a temperature gradient. The temperatures T° 2 and 
are found by application of the boundary conditions on 
the heat currents, which are proportional to d x T^ pm ^ s ^ 
with Tj = = T% if the heat fluxes have symmetric 
boundary conditions. 

Recall that T = To + arc for an isolated system under 
an applied temperature gradient. 



Using Eq. (29), substitution of Eq. (34 1 into Eq. (27) 



(which applies only to the inhomogeneous parts of 
T( s ,p,m)) gives, with no new parameters, 



Eq. (2). 



The coupling factor (R mp R ps < 1) between these 
modes further splits the two solutions; for R mp R ps 7^ 0, 
the (shorter) characteristic length A 2 = l/<?2 decreases 
and the (longer) length Ai = 1/qi increases. For three 
values of q m p/qps > 1 ; figure n3] shows the characteristic 
lengths Ai and A 2 , normalized by the pure mode phonon- 
magnon relaxation length (X mp = l/q mp ), versus the cou- 
pling factor R mp R ps . For q ps > q mp the plots are the 
same when Xi and A 2 are normalized by q ps rather than 

qmp • 



2 

£ 

7=1 



ips 



q 2 



d F K p 

x [T 7 a sinh fax) + T 7 b cosh fax)] 



T m = T + ax + 

7=1 



Imp 



\ Hp ) 



(35) 



x [T° sinh (g 7 x) + T 7 b cosh (g 7 x)] . (36) 



Substituting Eqs. (34), (35) and (K36J into Eq. ( 23 ) in turn 



gives the heat current in each subsystem: 



37 



37 



hip CO, tvp 

9 7 [T 7 Q cosh fax) + T 7 b sinh fax)] , 

7=1 
2 

n s a — k s g 7 [T" cosh (g 7 x) + T 7 sinh fax)] 



7=1 



^ <?7 [ r 7 Q cosh (97 1 ') + T 7 sinh («7 a; )] 

7=1 



'ips 



q 2 

imp 



d s ft s +df ftp \ 2 

<z 7 



a 2 

^rnp 



(37) 
(38) 

(39) 



r 



The total heat flux in the ferromagnet j J F = j7 + is 



37 



(k p + K m ) a - 



f i7 



7=1 



(k„, + K p ) g 7 ) 



a 2 

^imp 



I K, m +K p \ 



x [T 7 cosh (g 7 x) + sinh (g 7 x)] . (40) 



The boundary conditions on j a 



at a; 



-L/2 and 



a; = L/2 give a, Tf >2 and T^. 

Because heat flux is continuous, the total heat flux 
(integrated over all subsystems) due to each surface mode 
must be zero. This condition is satisfied by Eqs. (37), 
(38), and (39) on substitution from Eqs. (29) and (31). 



There are five unknowns in Eqs. (37), (38), and (39) 
(a, Tf , T$ , Ti, and T|), and seemingly six boundary 
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FIG. 3. The effect of mode coupling on the characteristic 
lengths associated with thermal equilibration in the spin- 
Seebeck system. The two characteristic lengths Ai and A2, 
normalized here by X mp = q mp , are shown as functions of 
the coupling factor RmpRps < 1, for: (a) q mp = q ps , which 
corresponds to equivalent pure mode lengths A mp = A ps ; 

(b) q m p = 3q ps , which corresponds to X pa — 3A mp ; and 

(c) q mp = 10g ps , which corresponds to X ps = 10A mp . For 
q ps > (Imp, the plots are the same when Ai and A2 are nor- 
malized by X ps rather than A mp . By definition, RmpRps < 1- 



conditions (for each of the three fluxes, one at x = —L/2 
and one at x = L/2). However, because the total en- 
ergy flux is conserved (i.e., no losses at the top of the 
ferromagnet dp or at the bottom of the substrate — d s in 



Fig. [TJ, there are only five independent conditions. 

For comparison to the theory of Ref. El we now con- 
sider the bulk system if the heaters contact the sam- 
ple and there is no substrate (so that g| 



?! 







lis)- 



1 2 m P and 



Then j e x F — > — (n p + n m )a, which repro- 
duces the homogeneous result of Ref. [6] (where Q = j e x F ) , 
and satisfies the condition of zero total heat flux due to 
the surface mode. If the heaters directly transfer en- 
ergy only to and from phonons (so that heat flow in the 
magnon system vanishes at x = L/2 and x = —L/2), 
then T£ — > K m a/[q mp K p cosh.(q mp L/2)] and T| — > 0, 
which reproduces the inhomogeneous solution of Ref. [HI 
As noted above, because T-f 2 are associated with a term 
proportional to sinh {q\^x) m the heat flux, = = T\ 
for symmetric boundary conditions on the heat fluxes 
(i.e., the same heat current is injected into each system 
at the "hot" side as is withdrawn from each system at 
the "cold" side). 

The above omits any consideration of how heat flows 
across the sample-substrate interface, which we now ad- 
dress. 



IV. HEAT FLOW IN 2D 

We have so far omitted any consideration of how heat 
flows across the sample-substrate interface, which is now 
addressed. We now consider heat flux along z, to explic- 
itly permit heat transfer between the substrate and the 
sample. We first detail the analytic theory, then present 
its numerical solution. 



A. Analytic Results 

To completely describe the z-dependence of the 
temperatures and heat fluxes in the system, the z- 
dependence of the heat flux input by the heater at 
x = —L/2 must be considered. In principle, it may have 
any functional form, and therefore properly requires a 
Fourier series in sin (fcz) and cos (kz) that includes an 
infinite number of lengths fc _1 associated with the z- 
direction. However, if the thickness (along z) of the sub- 
strate is much smaller than its length (along x), then k~ l 
should be very small compared to Ai 2 



q^l ofEq. (|31| 



The contributions from this z-dependence should decay 
along x over a distance on the order of the non-uniformity 
along z, and therefore we do not explicitly include them 
in the analytic theory. The cost of neglecting these high 
k values is that we cannot specify a heat input with a 
complicated variation along the thickness. 

We thus generalize equations ( 34 )-( 36 1 to take the form 



N 

+ E [ T (s, P ,m) n (*) sinh (q n x) + T b {s p m)n (z) cosh (q n x)] . 

n=l 

(41) 
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Note that we permit there to be N surface modes; for 
heat flow along only x, the one-dimensional heat equa- 
tions guarantee that N — 2, but the two-dimensional 
equations are nonlinear so that any N is allowed. 

The forms of Tf spm) Jz) and T^^Jz) are deter- 
mined by the conditions on the heat flux. We take sym- 
metric boundary conditions on heat flux along x, which 



give T, fe 



p,m) 



0. Then, substitution of Eq. (411 into 



Eq. ( 23 ) gives the heat fluxes along x and z to be 



3x 



N 



-(a,p,m) 



«( S ,p,m) 9» T (lp,m)„ ( z ) cosh (inX) , (42) 

71=1 

N 

K( s ,p,rn) d * T (s,P,m) n ( Z ) Sinh (inX). (43) 



This section finds the functional forms of T? B m x (z) and 
shows their amplitudes for example material parameters. 
It also discusses the bulk and boundary conditions that 
permit determination of their amplitudes, with the de- 
tails of these conditions given by Appendix [X] 

On properly treating the heat transfer between sample 
phonons and substrate phonons as z-directional currents, 
and heat transfer between sample magnons and sample 
phonons as a source/sink as for the ID case, employing 
Eqs. p3l and ((17) gives 



dX = 0, 



Cn 



;< !/ / - ; - /' — [T-m Tp)) 



(44) 
(45) 

(46) 



Equations (44|-(46) are identical to Eqs. (24)-(26l, but 



with phonon-substrate heat transfer in the form of fluxes 
rather than sources. These equations give 



^o m — To — Tq, 



(47) 



but they do not explicitly impose any conditions on T 0s 
or a s . For steady-state flow, however, we must take 



To. — T , 



(48) 



This relation guarantees that for any two of K( SiP , m ) to go 
continuously to zero, we recover the expected j% = —na. 



We now find T, a „ s (z) by substituting Eq. (411 into 



" (s,p,m) r . 



Eq. (|44|) and the decoupled forms of Eqs. (451 and (46) 



Substitution of Eq. ( 41 ) into Eq. ( 44 ) gives 
d 2 z T*(z) = -q 2 n Tl(z), 



(49) 



so that T" (z) is sinusoidal: 

Tl (z) = A£> cos (q n z) + A& sin (q n z). (50) 

Here, A^J and A s 2 } are constants determined by condi- 
tions on heat flux (see Appendix [A]) . 



Decoupled equations for T p and T m , and thus for 



T" (z) and (z), are found by combination of Eqs. ( 45 ) 



and (|46|). Addition and subtraction gives 

0, 



Kp&i T p -I- K m 3j T n — 2 (T m Tp 



Combination of Eqs. (|5T|) and (|52|) gives 

0, 
= 0. 



dld)Tp 



2 A^ry-? 

1mv a i 1 P 



d 2 d 2 T -a 2 d 2 T 



(51) 
(52) 

(53) 
(54) 



where we have employed Eq. ( 29 ) . Use of Eq. ( 41 ) in 



Eqs. (53 1 and (54) gives, for each mode n, 

^^(p,m)„ ( z ) + QnT(p,m) n ( z ) + ^Qn^zT{p,m) n ( Z ) 
~QmpdzT( l p,m) n ( Z ) ~ ?mp9n^~(p,ra)n (%) = 0. (55) 



The solution of Eq. (|55|) is 



T(p, m ) n ( z ) 



AY> . eVrt 



+ A 



(2) 
(P,n 



Here, A 



+ A fp]m) n COS + A fp]m) n Sm (inZ). (56) 

(1,2,3,4) 



are constants determined by conditions 



on heat flux (see Appendix \M 



Due to the mode splitting discussed in Sec. Ill the ID 



inverse lengths straddle q mp , that is, q^ 1 ^ > q mp > q^ . 
Therefore, for coupling R m pRps 7^ 0, the exponential 
terms in Eq. (56) are, in fact, oscillating terms for each 

mode that has q n > q^ D ^ ■ 



B. Bulk and Boundary Conditions 

Although T , a, A£' 2) , and A^Af are 2 + 10 A un- 
knowns associated with the temperatures and heat fluxes, 
they are not free parameters. As shown in Appendix \K\ 
bulk energy conservation gives AN conditions; energy 
conservation at the boundaries z — ~d s and z = dp, 
where we assume no heat loss to the vacuum, gives 3A 
conditions; there arc 2A conditions on heat flux at the 
interface z — 0; and there are 2 + N conditions on tem- 
perature and heat flux near the boundaries x = ±L/2. 
With these conditions, the present theory has no fitting 
parameters. 

Specifically, the 3A^ boundary conditions at z = — d s 
and z — dp are given by 



f z m (x,z = dp) = 0, 
j%»(x,z = dp) = 0, 
f z *{x,z = -d s ) = 0. 



(57) 
(58) 
(59) 



As discussed in Refs.|71IHl andHH heat currents are driven 
across an interface by the temperature difference across 
the interface, so that 

f/ {x, 2 = 0) = -h K [T p {x, z = 0)-T s (x,z = 0)] , (60) 



which gives N conditions. At the interface we take heat 
to be transferred only between substrate and sample 
phonon systems, so that 



f z "(x,z = 0) = f/(x,z = 0), 



(61) 



further below. For the numerical solution, we therefore 
estimate X mp = 2 mm from the observed voltage decay 
length in Fig. 2 of Ref. HJ We now present the results of 
this method, calculated using Mathematica v. 8.0. 



or equivalently 



f™(x,z = Q) = 0, 



(62) 



giving another N conditions. One imposes any two of 
Eqs. (60), (61), and (62), with the third being implicitly 



guaranteed by the energy conservation in the equations 
of motion. 

Only the remaining conditions, associated with the 
boundaries x = —L/2 and x — L/2, can be varied: 
the average temperature Tq, the temperature gradient 
a, and one condition per mode, associated with the rel- 
ative amounts of heat carried by each subsystem at a 
given short distance from the heater. All of these 2 + N 
conditions are set by experiment, the first two of which 
are, respectively, proportional to the sum and difference 
of the heater and heat sink temperatures. The other N 
conditions are related to the relative amounts of heat flux 
carried by each subsystem near the heater or heat sink. 
These conditions are non-obvious, but Appendix [X] ar- 
gues that they may be approximated by assuming that 
near the heater the heat flux carried along x by the sub- 
strate phonons dominates that carried by either the sam- 
ple phonons or sample magnons. 



TABLE I. Parameters used in numerical calculations, results 
of which are shown in Fig. [6] '"'Taken from Fig. 3 of Ref. [TO] 
' 6 'To our knowledge, this has not been measured, so we make 
an order of magnitude estimation. ^'Value unknown; K m /K p 
is likely to be lower at high temperature. ^'Estimate from 
Fig. 2 of Ref. [2] for the decay length of the observed spin- 
Seebeck voltage signal. '^Estimate for Rh:Fe on AI2O3 from 
Fig. 34 of Ref. U 



Parameter 



ftm I ftp 



d 3 
h K 



Value 



500 
100 
1/10 
1 x 10~ 7 
5 x 10~ 4 
5 x 10 2 
1 x 10 7 
15.5xl0" 3 



Units 



W/m-K 
W/m-K 

m 
111 

1T1- 1 

W/m 2 -K 
111 



Ref. 



(b) 
(c) 

M 
M 

Ef e) 

m 



Numerical Solution 



One can not assume that the inverse lengths for ID 



heat flow, given by Eq. (31) and now called q[ 1D ^ and 



q^ D \ are equivalent to the inverse lengths associated 
with 2D flow. Indeed, numerical solution with either 

of q[ 1D ^ or g^ 10 "* can be shown to be inconsistent with 
energy conservation. Since the 2D heat flow equations 
are nonlinear, analytic solution is not possible in general. 
However, an iterative approach can be used to find consis- 
tent values for q: solve the appropriate boundary condi- 
tions for the mode amplitude coefficients (i.e., the coeffi- 



in Eqs. (50) and (56)) using qi n i t = q[ 1D ^ 



cients A, 

{s.p,m) 

or (frnit = g , 2 1D ''; using these values for the coefficients, 
find the q new that guarantees energy conservation; begin 
the loop again using an appropriately chosen q' init in be- 
tween gi n j t and g new - One must iterate until q aew and 
ffinit converge P 1 

For our numerical calculations, we use the material 
parameters given in Table [TJ Note that Ref. [5] estimates 
A mp to be at least an order of magnitude too small to be 
the unusually large decay length of the observed voltage 
difference AV y , and the present theory does not explain 
such a large discrepancy, because as shown in Fig. |3j we 
do not predict mode coupling to amplify the larger length 
by a full order of magnitude. This matter is discussed 



Following Table [TJ Eq. ((MJ gives 
q[ 1B) = 476.73 m" 1 , q[ 



(1D) = 476.73 m -1 , oi 1D) = 1.0000 x 10 6 m -1 . 



(63) 



Using these as trial values for the numerical solution of 
2D heat flow boundary conditions, we find 2D inverse 
lengths consistent with energy conservation to be 



(2D) 



476.73 m" 



Although q 1 ^ ' and qX"' match to one part in 10 s (not 

shown to this precision above), only q^ 2 ^ satisfies energy 
conservation. 

The subsystem contributions to heat flow along z and 
along x for the two modes associated with g 1 2D ' ) and q^ D ^ 
are respectively shown in Figs. [4] and [5] Fig. [5| explains 
the significant difference between q£ and q^ D ^ ; the 
ID solutions q[ 1D ^ and q^^ should apply for heat flux 
along a; uniform in z. This holds for the q^ D ^ mode in 
5a, whereas the q^ D ^ displays significant curvature 



<k 



(2D) 



1.0015 x 10 b m~ 



(64) 



,(2D) 



D. Infinite Number of Inverse Lengths 



Other consistent solutions > q^ D ^ > q^ 2 ^ can be 
found numerically. We are here searching for the normal 
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- Kp d z T^(z) 



d z T° ( Z )x2xlO" 



FIG. 4. The phonon and magnon heat fluxes (in arbi- 
trary units) along z, for a given x, as a function of z, i.e., 
~> t (s,p,m)d z Tf l s p ^^(z), in the thermal equilibration modes 
with the two largest characteristic lengths. The substrate oc- 
cupies 2 < and the sample, with thickness magnified by 10 3 , 
occupies z > 0. In the sample the magnon heat flux is nearly 
parabolic and the phonon heat flux is nearly linear. In (a), 
where n = 1, the heat flux in the substrate is nearly linear. 
In (b), where n = 2, the heat flux in the substrate has many 
oscillations because A^ 20 ' >C d s . For both modes the sample 
is too thin for magnons to build up significant heat flux along 
z; in both (a) and (b) the magnon heat fluxes are magnified 



by 2 x 10 



10 





K P^P1 ( 


-d a 




dp 




-K m T^(z) x 10" 3 





b) 



s T«{z) 



d% 



FIG. 5. The phonon and magnon heat fluxes (in arbi- 
trary units) along x, for a given x, as a function of z, i.e., 
— K (s,p,m)Tu p m ) (z), in the thermal equilibration modes with 
the two largest characteristic lengths. The substrate occupies 
z < and the sample, with thickness magnified by 10 3 , oc- 
cupies z > 0. In (a), where n = 1, the magnon heat flux is 
multiplied by 10 -3 . For the parameters of Table [i] (a) shows 
that along x the heat flow for n = 1 is carried by all three 
subsystems, with magnon heat flux opposing sample and sub- 
strate phonon heat flow, and (b) shows that along x the heat 
flux for n = 2 is carried mostly by the phonon subsystems, 
which oppose one another at the interface. In (b), where 
n = 2, the heat flux in the substrate has many oscillations 
because A^ 21 "^ >C d s . Although it is not obvious at this scale, 
each heat flux has some curvature. 



modes associated with heat flow with the largest decay 
lengths, the larger q (and therefore smaller A) solutions 
are irrelevant to the current discussion. We do, however, 
discuss the nature of these solutions. 

Figure [7] shows the magnitude of the seven smallest 
wavevectors (except qi) versus the number of the solu- 
tion n (numbered by magnitude with q^+l > Qn^)- As 
n grows, the difference Sq between the inverse lengths of 
successive modes approaches either n/d s or 7r/(d s +cIf)', 
since d s 3> dp , it is difficult to distinguish which is the 
limiting quantity. Thus, the higher solutions are associ- 
ated with the geometry of the system. We do not discuss 
them further. 



Note that this numerical method, which searches for 
consistent values of q by using trial values, might not 
obtain all solutions, no matter how exhaustive the list of 
trial values. However, any missed modes are expected to 
have large q and small A, and thus are irrelevant to the 
current discussion. 



V. ON THE MEASURED EXPONENTIAL 
LENGTH 

For the calculated maximum X mp of Ref . [5J the present 
theory cannot account for the anomalously large length 
(on the order of 1 mm) observed in the spin-Seebeck ex- 
periments. On one hand, for the sample-substrate length 
Xp S to be on the order of 1 mm, with k s « k p ~ 10 2 W/m- 
K, d s ~ 100 nm, and dp ~ 10 nm, Eq. (29) gives an 
abnormally small thermal boundary conductance Hk ~ 
1 W/m 2 -K. Although Hk is not known for the particu- 
lar combinations of materials used in Refs. Fig- 34 
of Ref. E gives h K 10 7 W/m 2 -K (for Rh:Fe on A1 2 3 
at T — 50 K). We do not expect that thermal matching 
between substrate and sample in the spin-Seebeck exper- 
iments to be considerably worse. On the other hand, 
for the magnon-phonon length X mp to be on the order 
of 1 mm, the mode coupling term given by R mp R ps 



in Eq. (31) would have to account for a large increase 



of X mp (at least three-fold in the case of PermalloyP) 
Because spin-See beck experiments are carried out near 
room temperature^' or at T > 40 Kpl it is unlikely that 
the magnons carry a significant amount of the heat flux 



10 




FIG. 6. The relative magnitudes of phonon and magnon heat flux along z as a function of x and 2, i.e., jj s - p - m > hi arbitrary 
units. The substrate (only part of which is pictured) is at z < and the sample is at z > 0. The sample magnon heat flux is 
magnified here by the factor 3 x 10 11 ; for the parameter values of TablejlJ the sample is too thin for magnons to build up much 
heat flux along z. The profile of each subsystem's heat flux along z varies as sin (q n x). 



1.02 x 10 6 

1.00 x 10 6 
0.98 x 10 6 



77 



FIG. 7. The inverse lengths q n for n = 2 to n — 7. The 
inverse length q\, which is not shown, is three orders of mag- 
nitude smaller than 52- The difference 8q between the inverse 
lengths of successive modes quickly approaches a value near 
n/ds ~ ir/(d s suggesting that the additional modes are 

associated with the physical geometry of the system. 



in the ferromagnet, i.e., it is likely that K m <C k p . Since 
the mode coupling term R mp is proportional to K m /K p , 
mode coupling is likely a weak effect. 

However, phonon-magnon drag, as proposed in Refs. [HI 
and I10[ or some other mechanism may explain a much 
longer X mp than previously calculated. Hence, we have 
taken X mp to be larger than estimated by Ref. [5] (see 
Table |l]). The results above show that, for such a large 
X mp , m the spin-Seebeck system we expect a thermal 
gradient along z that varies as sinh (x/X), resembling the 
AV y measured by Ref. [5] (see its Fig. 2). 



VI. RELATING LONGITUDINAL THERMAL 
GRADIENTS TO TRANSVERSE VOLTAGE 
DIFFERENCES 

The relation between the applied longitudinal tem- 
perature gradient and the transverse voltage difference 
is complicated, and worth discussing. So far we have 
shown that the applied longitudinal temperature gradi- 
ent leads to a transverse (along z) temperature gradient 
in the sample - the first of th e three steps in Eq. ([!]), 



AT X EqmL > d z T. In Sec. VI A we show how to go from 
this tranverse temperature gradient to the accompany- 
ing transverse gradients of the magnetoelectrochemical 
potentials - the second of the three steps in Eq. ([I]), 

d z T — ^> d z prf,± - which are defined below. Finally, in 
Sec. |VIB| we show how to go from these transverse gradi- 
ents (along z) of the magnetoelectrochemical potentials, 
via the up- and down- spin Hall conductivities, to the 
measured transverse (along y) voltage difference AV y - 

the third of the three steps in Eq. ([T]), d z pf,i ISHE > AV y . 

We do not consider the use of platinum bars, which 
introduces a very complex geometry and is beyond the 
scope of the present work (and, as noted above, the effect 
has been observed with point contacts). 



A. On Magnetoelectrochemical Potential, 
Temperature, and Spin Current 

By irreversible thermodynamics, the total spin flux 
(defined below as the difference of the number fluxes 
of up- and down-spin carriers), is driven both by gra- 
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dients of temperature and of magnetoelectrochemical 
potentials The magnetoelectrochemical potentials^ 
are defined bjffSEl 



M. 



(65) 



Here, m and /ij, are the chemical potentials of up- and 
down-spin electrons, e is the electron charge, <fi is electri- 
cal potential, g is the electron ^-factor, fig is the Bohr 
magneton, H* is the effective magnetic field, and M is 
the direction of magnetization. The field H* is the dif- 
ference between external magnetic fields and the internal 
fields, including the exchange and dipole contributions, 
and is defined so that H* = in equilibrium. A more 
detailed discussion of H* is given in Ref. [1IJ] 

The up- and down- spin fluxes are primarily driven by 
the r especti ve gradients and p,±, but each has cross- 
^ erm gili6ii9| asS ociated with the other potential, as well 
as with the temperature. We thus write 

jj = -ifs^Tm - -J^iMt _ L U 9 iPl' ( 66 ) 



3i = -L l£ diT m 



<7| 

-2^1- 



(67) 



Here, cr-f and are the respective bulk conductivi- 
ties of up- and down-spins (generally not equal in a 
ferromagnet), and L^ £ and L^ e (with units of m/K-s) 
are cross-term coefficients associating thermal gradients 
and individual spin-carrier currents (thus associated with 
both the electrical and spin currents). By an Onsager 
relation 2 ^ £^ = (with units of m/J-s) are cross-term 
coefficients associating up and down spin currents with 
down and up magnetoelectrochemical gradients. Typi- 
cally L-fj, = are taken to be small, so that the terms 
L^diPl and L lt d iP"t are negligible. 

To calculate jj and jf everywhere, we employ their 
boundary conditions (that they have zero normal com- 
ponent at each sample boundary, which assumes no sur- 
face scattering) and their bulk equations, given for steady 
state by 



^j\ 



S 



S 



It- 



(68) 



For charge conservation, the up- and down-spin source 
terms S^i and (which are proportional to (u+ ~ 
f L l)/ T sf^ where t s j is a characteristic spin-flip tirn d 16 ! 19 ^ 
are equal and opposite. Substitution from Eqs. (|66|) and 



(67) into Eq. (68) gives two equations for two unknowns, 



fi-f and /ij,. Because the temperatures are shown above to 
vary as sinh (x/A), then p,^ and JX^ also vary as sinh (x/A). 



B. On the Spin Hall Effect 

We now discuss how to go from d z \x^ and d z fii to the 
measured voltage difference along y, i.e., AV y . We work 
by analogy to the Hall effect, which occurs when an elec- 
tric flux J is driven through a conductor in the presence 



of a magnetic field B' that is perpendicular to the cur- 
rent . 

Consider a conductor of width w along y. Let the 
electric current be driven along z by an applied electric 
field E z , so that charge carriers have a velocity v z . With 
an applied magnetic field (B' x ,0,0), a Lorentz force then 
drives the charge carriers along y, so that charges of op- 
posite signs accumulate at the edges. The Lorentz-force- 
induced current is given by Jj — av z B' x . In the steady 
state, there is no flow along y, so an electric field E y 
develops to oppose the Lorentz-induced current along y. 
The total charge flux along y is given by 



J y = = a (Ey 



*B' X ) 



The so-called Hall field E y thus is given by 



E„ 



-v z B' 



J Z B' 



lie 



(69) 



(70) 



where we have used J = —nev, and n and — e are the re- 
spective concentration and the charge of the charge car- 
riers. The Hall voltage is AV y = E y w. 

Thus, the Hall effect relates an applied electric current 
to a measured transverse electrical potential difference. 
In contrast, the Spin Hall effect (SHE) relates an applied 
electric current to transverse differences in the magne- 
toelectrochemical potentials, and the inverse Spin Hall 
effect (ISHE) relates an applied spin current to a trans- 
verse difference in electrical potential (see, for example, 
Refs. [2TH26]). For the SHE and ISHE there are fluxes of 
charge carriers with both up- and down- spin. Instead 
of the action of Lorentz force in the Hall effect, for the 
SHE there are forces due to the spin-orbit interaction, 
whose effect enters via non-zero up- and down- spin Hall 
conductivities <t s h+ and <J S H t - (Thus the effect of the 
spin-orbit interaction is taken to be a perturbation.) In- 
stead of the electric field E y = —d y 4>, the spin-orbit force 
is associated with —d y p,^ and —dypi- We take the con- 
tributions to the number fluxes along y of the up- and 
down- spin carriers by this spin-orbit force^ to be given 
by 



Jy — O z fXf, Jy 



sHl = (J eH l 

i e 



dzfii, (71) 



The number fluxes along y of the up- and down- spin 
carriers can be written as 



Jy = - — dyVt 

4 = -v^ 



ffsff t 

e 



(72) 
(73) 



For no charge current along y, the sum jy+jy = 0. We 
also assume no bulk spin current along y, so jj — jj = 0. 



Thus, we take jj = and jj = 0, so that Eqs. |72j) and 



(73) give 



dyfPf 



(J sH i 
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The known sources d z a^ and d z u^ on the right-hand- 
sides (RHS) of Eq. (74) are uniform in y. 



With AV y = dyS(j>, integration of Eq. (|82| over y 

across the width of the sample then gives 



To write the magnetoelectrochemical potential in 
terms of the concentrations of up- and down-spins and 
the electric potential, we linearize the chemical poten- 
tials and the effective magnetic field term as 



<Vf,. 



dn 



1M 



<W,4 



SH* . M = {5n t - Sni) , 

(75) 



where 6 denotes deviations from equilibrium, jiQ is the 
permeability of free space, and \ is the magnetic suscep- 
tibility. Then Eq. ( 65 1 gives 



dn- 



t4 



■<)n^ - e8cj> ± 9 -^- (Sn t - 6 H ) . (76) 



With dfif i /0n+ i uniform in y, substitution of Eq. ( 76 ) 
into the left-hand-sides (LHS) of Eq. (74 1 gives 



^^d y Sn t - ed y 8(j) + 9 ^^ B (d y 5iVf - dySn^) 



^^dySni - ed y 5((> - 9 ^^ B {d y 8n^ - d y 5n^) 



(77) 



(78) 



With the RHS of Eqs. (JT7J) and (|78j> known, they are two 
equations for the three unknowns 8n^, 8n^ 1 and 6<f>. A 
third relation is provided by Gauss's Law: 



e e 



{Sn t + 6ni) . 



(79) 



where £q & n d £ are the permittivity of free space and the 



relative permittivity. Solving Eqs. ( 77 )-( 79 ) gives 8n^, 



8n^, and 8cf>, the last of which is related to the measured 

voltage by AV y — f™^ 2 dyScj). We now discuss the solu- 
tion. 

It is consistent to take 8ru = —6n±, i.e., local 
electroneutralityP^ (Equations (|77])-([79| then give that 
d v 6(j) and d y 8n^ are uniform in y.) Equations (77) and 



781) can then be solved for d y 8n^ and d y S(f>. Defining the 



dimcnsionlcss ratio 



dn^ dn^ 



du t dm 2gu Q u 2 B 



(80) 



dn^ dn 



-I 



we have 



d y 8n^ = R 



\dn^ 
2e 



dn x 



(T-|- 



-d z fj^ 



a sH i 



dzjJ-t 



1 + R u \ a sH 



2e 



(81) 
(82) 



1) -^d z n t - + 1) -^0,04. 



<7-|- 



°4 



(83) 



where we have employed the uniformity of d z p^ and 
d z fli along y. As discussed above, j~L^ and ^ vary as 
sinh(x/A), thus Eq. (83) predicts AV y ~ sinh(x/A). 



The present work shows that the relation between AV y 
and AT X is very complicated, and suggests that a direct 
relation AV y ~ SsAT x (see, e.g., Ref. [3J is correct, but 
may not be quantitatively useful. However, the present 
work does support such a qualitative analysis, where the 
applied thermal gradient along x leads, via the spin- 
Seebeck effect, to spin carrier fluxes along z, which in 
turn produce the measured voltage difference AV y along 
V- 

When surface scattering occurs, the present analysis 
would become much more complex; see Landauer and 
SwansorP2J for the ordinary Hall effect. 



VII. SUMMARY AND CONCLUSION 

The present work finds the detailed temperature pro- 
file for the spin-Seebeck system, including both sam- 
ple and substrate, when a temperature difference AT X 
is applied along x. For ID heat flow (only along x) 
we find that the temperature contains a part varying as 
sinh (x/A), for each of two characteristic lengths (A ps and 
A mp ), one of which may correspond to the observed de- 



cay length of AV y . Equations (31) and (29) show that 
quadrupling the thickness of both the sample and sub- 
strate should approximately double these lengths. Pol- 
ishing (roughening) the substrate before depositing the 
sample should increase (decrease) hx , and thus decrease 
(increase) X ps . If X ps corresponds to the observed expo- 
nential decay length, measurements on a series of sam- 
ples with increasingly rough sample/substrate interfaces 
should reveal this dependence. Further, changing the 
coupling factor between the modes (by changing K m / 'k p 
or d s K s /dpK p ) modifies both lengths - increasing either 
increases the larger length, which likely corresponds to 
the measured decay length of A V y . 

For 2D heat flow (along both x and z), we also find 
that the temperature and thermal gradients along z in 
the spin-Seebeck system vary as sinh (x/A), and find a 
complicated sinusoidal and exponential profile along z for 
the thermal gradients, with an infinite number of charac- 
teristic lengths, which we study numerically. The longest 
of these corresponds to the longest ID length. The sec- 
ond longest length is a geometry-modified version of the 
other ID length. Further lengths are largely due to the 
geometry. 

We show how the thermal gradient along x leads to 
the measured AV y . The thermal gradient along x leads 
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to a thermal gradient along z, which then drives up- and 
down- spin currents along z (the spin-Seebeck effect), and 
is accompanied by gradients along z of the magnetoelec- 
trochemical potentials. These magnetoelectrochemical 
potential gradients along z then produce the measured 
AV y , via the inverse Spin Hall effect (due to a nonzero 
spin-orbit interaction that leads to spin-Hall conductivi- 
ties). 



VIII. ACKNOWLEDGEMENTS 



We would like to acknowledge V. Pokrovsky for valu- 
able conversations, and the support of the Department 
of Energy through grant DE-FG02-06ER46278. 



* wsaslow@tamu.edu 

1 M. Johnson and R. H. Silsbee, Phys. Rev. B, 35, 4959 
(1987). 

2 C. M. Jaworski, J. Yang, S. Mack, D. D. Awschalom, J. P. 
Heremans, and R. C. Myers, Nature Materials, 9, 898 

(2010) . 

3 K. Uchida, S. Takahashi, K. Harii, J. Ieda, W. Koshibae, 
K. Ando, S. Maekawa, and E. Saitoh, Nature, 455, 778 
(2008). 

4 K. Uchida, J. Xiao, H. Adachi, J. Ohe, S. Takahashi, 
J. Ieda, Y. K. T. Ota, H. Umezawa, H. Kawai, G. E. W. 
Bauer, S. Maekawa, and E. Saitoh, Nature Materials, 9, 
894 (2010). 

5 J. Xiao, G. E. W. Bauer, K. Uchida, E. Saitoh, and 
S. Maekawa, Phys. Rev. B, 81, 214418 (2010). 

6 D. J. Sanders and D. Walton, Phys. Rev. B, 15, 1489 
(1977). 

7 G. L. Pollack, Rev. Mod. Phys., 41, 48 (1969). 

8 E. T. Swartz and R. O. Pohl, Rev. Mod. Phys., 61, 605 
(1989). 

9 H. Adachi, K. Uchida, E. Saitoh, J. Ohe, S. Takahashi, 
and S. Maekawa, Appl. Phys. Lett., 97, 252506 (2010). 

10 C. M. Jaworski, J. Yang, S. Mack, D. D. Awschalom, R. C. 
Myers, and J. P. Heremans, Phys. Rev. Lett., 106, 186601 

(2011) . 

11 J. M. Ziman, Electrons & Phonons (Oxford University 
Press, London, 1960). 

12 Y. G. Gurevich and O. L. Mashkevich, Physics Reports, 
181, 327 (1989). 

13 1. 1. Hanna and E. H. Sondheimer, Proceedings of the Royal 
Society of London, Series A, Mathematical and Physical 
Sciences, 239, 247 (1957). 

14 F. Keffer, in Encyclopedia of Physics: Ferromagnctism, 
edited by S. Fliigge and H. P. J. Wijn (Springer, Berlin, 
Germany, 1966) p. 1. 

15 H. B. Callen, Thermodynamics (Wiley, New York, 1960). 

16 M. R. Sears and W. M. Saslow, unpublished. 

17 Care must be taken in determining a new initial value for 
the next iteration. For gi n i t far from a consistent value (that 
is, a value that satisfies energy conservation), then qimt 
and q ncw will differ significantly. Naively choosing g;' nit = 
9init + !(<2ncw — 5init) can result in a non-converging series. 
Hence, we include the factor 1/C in place of 1/2 to define 
ql n it', depending on the initial choice of qmit, convergence 
can require C ~ 10° or greater. 

18 ReferenceElappears to define the "spin potential" to be the 
difference of the up and down spin chemical potentials, and 
thus does not include the contribution from the magnetic 
disequilibrium associated with spin accumulation; this is 
different from the magnetoelectrochemical potential. 

19 W. M. Saslow, Phys. Rev. B, 76, f 84434 (2007). 



L. Onsager, Phys. Rev., 37, 405 (1931). 

M. I. Dyakonov and V. I. Perel, Phys. Lett, A35, 459 

(1971). 

J. E. Hirsch, Phys. Rev. Lett., 83, 1834 (1999). 

S. Zhang, Phys. Rev. Lett., 85, 393 (2000). 

J. Sinova, D. Culcer, Q. Niu, N. A. Sinitsyn, T. Jungwirth, 

and A. H. MacDonald, Phys. Rev. Lett., 92, 126603 (2004). 

E. Saitoh, M. Ueda, and H. Miyajima, Appl. Phys. Lett., 

88, 182509 (2006). 

E. M. Chudnovsky, Phys. Rev. Lett., 99, 206601 (2007). 
To eliminate the applied magnetic field B x , rather than the 
spin-orbit interaction, as the source of the deflection of the 
spin carriers leading to AV y , Ref. [4] shows that AV y = 
when Cu bars (with weak spin-orbit interaction) are used 
for detection, whereas AV y 7^ for Pt bars (with much 
stronger spin-orbit interaction). 

Because Ref. [3] takes the "spin potential" to be the differ- 
ence of up- and down-spin chemical potentials, it neglects 
the effective field H* and therefore by Eq. ( 75 1 assumes 



no spin accumulation. The present work, however, permits 
spin accumulation. 
29 R. Landauer and J. Swanson, Phys. Rev., 91, 555 (1953). 



Appendix A: Bulk and Boundary Conditions for 
Heat Flow along x and z 



With Eqs. (47) and (48) relating the linear terms in 
temperature, there are 2 + 107V unknowns in Eqs. (42), 



u 4(1,2) 
given by ' 



(|43|, (|50p and (|56|): one T , one a, and lOiV amplitudes 

(1,2,3,4) 



J 4.pi' 2 ' 3 ' 4 \ and Ami?' '' 1 ' . This section 



details the bulk and boundary conditions on heat flux 
that give these unknowns. 



1. Bulk Conditions 

By matching coefficients of like terms, substitution of 
Eqs. (Jill and (pel) into Eq. (1521) gives 



A (i) = 

m n Pn ' 

A (S) =A (3) 

Win Pn ' 



A< 2 > = 



1L 4(2) 

Pn 



m n Pn 



(Af) 



Since each of the above relations is a single condition 
for each mode n = 1, 2, ...N, then Eq. (Al) gives 4iV 
conditions. 
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2. Boundary Conditions 

a. Boundary Conditions on Heat Flux along z 

There are a further 5N conditions given by the bound- 
ary conditions on the heat flux along z for the various 
subsystems at z — — d s , z — 0, and z — dp. They are 
given above as Eqs. (57)-(59), and any two of Eqs. (60)- 



( 62 ) with the third imp 



icitly guaranteed by energy con- 



servation. 



b. Boundary Conditions on Heat Flux along x 

Two further conditions that constrain the homoge- 
neous temperature coefficients, To and a, come from the 
temperatures of the heater and the heat sink. The re- 
maining conditions on heat flux along x are not obvious. 

With the heater and heat sink each in contact only 
with the substrate, we take the boundary conditions in 
the x-direction on each energy flux j e x are symmetric (we 
employ this above in taking Th pm Jz) — 0). This pre- 
cludes permitting the heat flux input by the heater to 
have a different profile in z than the heat flux output 
to the heat sink. However, as stated above, we are only 
treating the region far enough away from the heaters that 



the details of heat flux entering and leaving at x = ±L/2 
are irrelevant. Only a full solution with an infinite sum 
over inverse lengths q n can treat the specifics of the in- 
terfacial input, and it is beyond the scope of this work to 
solve for infinite inverse lengths. Thus, we can not apply 
boundary conditions precisely at x = ±L/2. 

We make the following approximation: at x = ±L/2=F 
is, where is is just far enough away from the heater/heat 
sink that the details of the input heat flux are irrelevant, 
we take d x T p = and d x T m = 0. We take the heaters to 
be in contact only with the substrate, and assume that a 
significant amount of heat does not seep into the sample 
over the distance is- Explicitly, 



^a;T m [X 
&x ( X - 



-L/2- 
-L/2- 



■ls) 
is) 



= 0, 
0. 



(A2) 
(A3) 



Recall that we take heat flux (and therefore d x T) to 
be symmetric about x — 0, so that the conditions at 
x = +L/2 — is are not independent. Although it is not 
obvious, Eqs. (A2) and (A3) give N conditions, which 



relate the amplitudes of each of the N surface modes to 
the others. 

Thus, for the 2 + lOA^ unknowns in the substrate 
phonon, sample phonon, and sample magnon temper- 
atures associated with heat flow along both x and z, 
Eqs. ||57|)-([62]) and (|ATj)-(lA3l give 2 + 10 AT conditions, 



and there are no free parameters. 



